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Introduction 

Hawking radiation is one of the most spectacular and unexpected predictions of mod- 
ern theoretical physics. In 1974 Hawking showed by combining General Relativity and 
Quantum Mechanics, that stationary (even static) black holes are expected to emit thermal 
radiation at a temperature proportional to the surface gravity of their horizon. This had 
a profound impact on theoretical physics leading to a beautiful synthesis between gravity 
and thermodynamics, cathalized by quantum mechanics. Unfortunately, the astrophysical 
relevance of Hawking radiation is negligible: solar mass black holes should radiate at a tem- 
perature of the order of 10~^ K, much below the cosmic microwave background. There is no 
hope to identify such a tiny signal in the sky. Nevertheless there is no doubt that Hawking's 
result remains a milestone in the quest of unifying gravity and quantum mechanics. 

The basic feature of Hawking radiation is that it is the result of a pairs production 
process which proceeds by the conversion of vacuum fluctuations into on shell particles. The 
possibility of having this pairs production in a stationary or even static setting lies in the 
existence of negative frequency modes. These indeed exist inside a BH horizon. They are 
trapped modes which are associated with the so called "partners" of the Hawking thermal 
quanta, that reach the asymptotic region far from the black hole ^ . 

In Fig. 1 there is a sketch of the light cones in a black hole (BH) spacetime and in 
Fig. 2 of the modes involved in the Hawking process. These are outgoing, i.e. in principle 
right-moving. However inside the BH horizon they are unable to escape and are forced by 
the strong gravitational field to move to the left towards the internal singularity. 

Note that the singularity plays no role at all in the Hawking process. It is just the 
future end-point of any physical trajectory inside the BH: it has no causal effect on the 
spacetime. One should also mention that the ingoing (left moving) modes are basically 
unaffected by the presence of the horizon, and participate only marginally to the process 
through backscattering. 

Finally, Hawking radiation does not depend on the details of the formation of the BH 
(stellar collapse for instance). These determine only an initial transient emission which 
rapidly (typically in an exponential way) decays in time, leaving a steady (ignoring backre- 
action effects) thermal radiation. This is Hawking radiation, which depends only (as required 
by the "no-hair" theorem) on mass, charge and angular momentum of the BH through the 
surface gravity 

The most serious drawback in Hawking's result is the so called "transplanckian" problem 
[sl. Because of the exponential redshift suffered by particles modes constituting Hawking 
radiation in their journey from the horizon region to infinity, only very high energy modes 
(whose frequency is much larger that the Planck one) created near the horizon manage to 
reach the asymptotic region. This fact makes people rather uneasy since it requires an 
extrapolation of our knowledge of physics to an unkown land characterized by scales less 
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BLACK HOLE 




Figure 1: Light-cones in a black hole spacetime 




HORIZON 

Figure 2: Modes involved in the Hawking process 

than 10~^^ cm, where quantum gravity effects might in principle completely alter the picture. 
We shall come back to this point later. 

We move now to white holes (WH). These are a kind of time reversal of BHs. In Fig. 3 
the propagation of light (light cones) in a WH is depicted. One sees that the horizon once 
again acts as a semipermeable membrane. For a BH (Fig. 1) it is permeable only from 
outside to inside. For a WH (Fig. 3) the horizon is permeable only from inside to outside. 

While a BH swallows everything which falls inside the horizon, a WH expels everything 
that lies inside the horizon, the ultimate origin of this matter being the singularity. Note 
that nothing can penetrate the WH horizon from outside. Ingoing physical trajectories pile 
up along the horizon which is an infinite, this time, blueshift null surface. This causes WHs 
to be unstable and therefore of limited interest in gravity as compared to BHs 0] . 

Nevertheless one can infer that, beacuse inside the WH horizon there are negative energy 
modes, there should exist a sort of Hawking radiation also for them. The modes now involved 
are the ingoing (potentially left moving) ones, as can be seen in Fig. 4. Here we see that 
the transplanckian problem appears in the future since the relevant modes pile up along 
the horizon and become infinitely blueshifted. Unlike the BH case, the process seems to be 
highly sensible to initial conditions. 

A much more serious problem is the presence of the singularity which now lies in the past 
and causally affects the rest of the spacetime, making every physical prediction about WHs 
quite arbitrary. 

As we have already stressed, at the heart of the Hawking effect is the existence of negative 
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Figure 3: Light-cones in a white hole spacetime 

WHITE HOLE 





Figure 4: Hawking process in a white hole spacetime 



energy states which allows pair creation out of the vacuum even in a stationary setting. 
These peculiar states exist inside the horizons of both BHs and WHs. But this is not at all 
an exclusive feature of these exotic systems. It can be found in many other "more normal" 
systems of different nature, most notable in fluids whose motion becomes supersonic in some 
region jsl. 

One can envisage two flows configurations which are the fluid analogues of the gravita- 
tional BH and WH, see Fig. 5. It is easily realised that sound waves propagation in the two 
fluid configurations proposed has the same qualitative behaviour as light propagation in a 
BH and a WH respectively. One indeed sees that in Fig. 5a sound waves in the supersonic 
region are trapped and dragged to the left by the flow as it happens for light inside a BH 
(see Fig. 1). Fig. 5a represents a sonic BH and similarly Fig. 5b represents a sonic WH, 
in which sound waves in the supersonic region are forced to move towards the subsonic one. 
The place where the flows turn from subsonic to supersonic (or viceversa) plays the role 
of the horizon (sonic horizon). The existence, here, of a trapped region bounded by the 
horizon, unlike the gravitational case (singularity theorems 6|), does not imply the presence 
of a singularity, and this makes in particular the study of sonic WHs physically meaningful. 

Because there are negative energy states inside sonic horizons, one expects the presence 
of an analogue Hawking radiation (of phonons) in both sonic BH and WH. 

The "transplanckian problem" associated to Hawking radiation translates in the fluid 
context by the fact that the relevant sound modes have characteristic wavelengths smaller 
than the fluid constituents distance, which clearly makes no sense at all. While there is 
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(b) WHITE HOLE 
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Figure 5: Flow configurations analogue to a BH (a) and a WH (b). 



no hope to address properly the transplanckian problem in gravity (we do not have yet a 
consistent quantum gravity theory), in the fluid case one often has a complete and well 
tested microscopic quantum mechanical description of the system at the molecular and 
atomic scale. 

The very existence and related features of the analogue Hawking radiation can then 
be investigated at a fundamental level and the theoretical predictions hopefully tested in 
laboratory experiments. The most favourable experimental setting to this end is probably 
offered by Bose-Einstein condensates 01, where the generally huge difference between the 
intrinsic temperature and the Hawking one can be signiflcantly reduced. 

In BECs the dispersion relation ceases to be linear at high waves momenta and becomes 
"sup er luminal" . The behaviour of the relevant modes involved in Hawking radiation in a 
BEC signiflcantly differs from the one corresponding to the gravitational case and depicted 
in Fig. 2 (BH) and Fig. 4 (WH). Indeed for the BEC BHs and WHs conflgurations one 
flnds, respectively, the behaviours in Fig. 6a and Fig. 6b. 



(a) BLACK HOLE 



(b) WHITE HOLE i 
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Figure 6: Behavior of the Hawking quanta and the partner in a BH (a) and a WH (b) in the case 
of a supersonic dispersion relation. 



One clearly sees that there is no longer an inflnite blueshift (in the past for BH, in 
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the future for WH), the modes no longer pile up along the horizon and no exponentially 
large frequencies are involved. The "transplanckian problems" simply do not exist for BEC 
acoustic BH/WHs, the modified dispersion relation has completely removed it. 

The aim of this paper is to provide a detailed study of Hawking-like radiation in BEC 
undergoing a WH like supersonic flow jsl. To allow an analytical (rather than numerical) 
treatment, which leads to a clear insight in the physical process responsible for the emission, 
an idealised setting is considered in which a step-like discontinuity in the sound velocity 
separates a subsonic homogeneous region from a supersonic one [ol, [lo|. Similarity and 
differences between BH and WH emission will be outlined. 



Stationary stepwise homogeneous condensates 
A. General framework 



The theoretical setup we shall use is the same adopted in ref [10[ for discussing Bogoliubov 
phonon creation "a la Hawking" in a BEC undergoing a BH like flow. We shall outline here 
the basic ingredients. 

A system of interacting bosons, confined by an external potential Vext, is described in a 
second quantized formalism by a field operator ^(t, x) which annihilates an atom at t, x. 
^' satisfies the usual bosonic commutator rules 

[^(t, x), ¥{t, x')] = 6'^{x - x') . (1) 

The time evolution of \l/ is given in the diluite gas approximation by 

ihdt^ = (^-|^V' + Vext + 9^^^^ ^ , (2) 

where m is the mass of the atoms and g is the nonlinear atom-atom interaction constant. 
At low temperatures (T ~ lOOnK), most of the atoms condense in a common ground state 
which is described by a c number wave function \E'o that satisfies Gross-Pitaevski equation 

ihdt^o = ("I^V^ + Kxi + gn^ *o , (3) 

where n = |\&oP is the number density of the condensate. Linear perturbation around this 
classical macroscopic condensate can be parametrized by an operator (p as follows 

^ = *o(l + 0). (4) 
satisfies the Bogoliubov-de Gennes equation 

^h^t4> =-( + -^V 1 + mc2(0 + 0t) , (5) 



2m m \1/ 
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Figure 7: Acoustic white hole configuration made of two semiinfinite homogeneous condensates: 
one supersonic (x < 0) and one subsonic (x > 0). 



the speed of sound, a fundamental quantity in the 



where we have introduced c = W 
following discussion. 

The configuration we shall consider is very idealized, consisting of two semiinfinite one 
dimensional stationary homogeneous condensates attached along a step-like discontinuity. 
In one region, say the left (x < 0) the flow is supersonic, while in the right one (x > 0) the 
flow is subsonic. The discontinuity is located at x = 0. See Fig. 7. 

The condensate is supposed to flow from left to right at an everywhere constant velocity 
Vo- This configuration, as mentioned in the Introduction, resembles a white hole, since sound 
waves (i.e. hydrodynamical excitations of the BEC) in the supersonic region are not able to 
propagate upstream, they are dragged by the flow towards the horizon (x = 0) and ejected 
in the subsonic region. Moreover, sound waves in the subsonic region are unable to penetrate 
the horizon which acts for them as a repulsive surface. 

Furthermore we assume that the condensate has everywhere (i.e. in both regions) a spatial 
uniform density n. One can however have different (constant) sound velocities in the two 
regions (say, q and Cr) by changing accordingly the interaction coupling g (remember that 
c = using the so called Feshbach resonance method 3]- If one changes the external 

potential from one region to the other as 



V 

9 n , 



(6) 



vLt + = const. = v;;, 

where /(r) refers to the left (right) region, one has that a simple plane- wave 

^o(a;,t) = v^e*'^"™"* (7) 

with ^ = vq and wo = ^ + {gn + Vext) is for all time a solution of the Gross-Pitaveski 
equation (E]) for every x. 

Concerning the (non hermitean) fluctuation operator (f), owing to the stationarity of our 
configuration, it can be decomposed in a positive and negative frequency part, namely 



(8) 
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and a] are phonons' annihilation and creation operators respectively. They satisfy bosonic 
commutation relations [aj,aj] = 6ij. 

Inserting the decomposition ([8]) into the Bogoliubov-de Gennes equations ([5]) and its her- 
mitean conjugate, we see that the modes satisfy the coupled equations 



i{dt + vod.^) + —dl - - 



(9) 



where ^ = h/mc is the so called healing length, which, as we shall see, sets the fundamental 
scale in our system, it is the analogous of the Planck length in gravity. 
As said above, because of the stationarity of the flow, the time dependence for both (pj and 
(fj can be taken as e~™*, where w is the conserved frequency of the modes as measured by 
a static (laboratory) observer. So w will label our modes, i.e. we shall work with modes at 
fixed w. The modes are normalized according to 



5„ 



hn 



(10) 



where the +(— ) sign refers to positive (negative) norm states. 

Piecewise homogeneity allows in each region, r or Z, the solutions for the modes to be given 
in the form of plane waves 



E(w)e 



-iwt+ikx 



(11) 

(12) 



where D{w) and E{w) are normalization factors to be determined using eq. (|T0|) . Inserting 
the plane wave solutions in the Bogoliubov- de Gennes equations we obtain 



w — vok) — 
-{w - Vok) - 



ice 

2 

ice 



E{uj) 



(13) 



This homogeneous linear system has nontrivial solution if the associated determinant van- 
ishes. This gives the dispersion relation 



{w - vokf = ^[e + 



ee 



(14) 



which can also be rewritten as 



w 



Vok = ^±{k) = ±cy e 



(15) 
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Q±{k) is the (nonconserved) frequency of the modes as measured in a frame comoving with 
the fluid. The ib sign in eq. f|T5|) labels the positive (+) and negative (— ) branches of the 
dispersion relation. The normalization condition fITU]) gives 

dk 



dw 



(16) 



which, with the help of eqs. (fT3|) . fixes the normalization coefficients as 

n( \ u-vok+'-f- 
D(uj) = 



Anhnc^k^ 



B,.) . - , . (17) 



^A7rhnc^k^\{uj-Vok) (f) 



Using these explicit forms of the normalization coefficients, one realizes that positive norm 
states belong to the positive branch of the dispersion relation, negative norm states to the 
negative branch. 

From eq. ( 1151) one can also see that to a w > negative norm solution of eq. (Q, one can 
associate a positive norm solution with w < (and k — j- —k). One can then work with 
only with w > modes, identifying w > negative norm states (if present) with negative 
frequency excitations. 

In the most familiar situation of a subsonic flow, for w > one has only positive norm states 
(negative norm states have w < 0). However, when the flow is supersonic we will see that 
for a given w; > there can be, in addition to positive norm modes, also negative norm ones 
which are physically interpreted as described above. 

Coming back to the dispersion relation, eq. (fT5|) . we see that for small k (i.e. A; ^ 1/.^) the 
dispersion relation is linear w — v^k = ±ck. This is the hydrodynamical regime {C, — )■ 0) 
with solutions ku = K = leading to Vg = ^ = vo±c for the group velocity of the 
two modes. In general, at fixed w, eq. ( !T5|) is fourth order in k, which admits four roots kw^ 
{i = 1, 2, 3, 4). So we can write the general solution of the Bogohubov-de Gennes eq. Q in 
each region as a linear combination of four plane waves with the above wave vector kw\ i.e. 

4 



0. = e-^^-'J2A^^D,{w)e^'^ , 
1=1 

4 

= e-™*^A;'")E.(i/;)e^'=''^ , (18) 

i=l 

where the Ai{w) are the amplitudes of the modes, not to be confused with the normalization 
coefficients Di{w) and Ei{w) which are fixed by the normalization condition and are given 
by eqs. f|T7|) with k = kw\ Note that the amplitudes Ai{w) are the same for both 0^ and 
ip^u as required by the equations of motion . 
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Figure 8: Dispersion relation in the subsonic (right) region. 

Once the roots of the dispersion relation are found in each region, one can construct the 
corresponding solution, eqs. f ITH]) . in the left and right region with the corresponding left 
and right amplitudes A[{w), Al{w). Note that the explicit form and features of the four 
roots will depend on the sub or supersonic character of the flow and so will be significantly 
different in the two regions. 

Now, the equations of motion, eqs. Q require and (f and their space derivatives to be 
continuous across the step discontinuity at x = 0, namely 

[0] = 0, [0'] = 0, [^] = 0, [ip'] = 0, (19) 

where [/(a^)] = lime_^o[/(2^ + ^) ~ /(^ ~ ^)] ^ prime means The four matching 
conditions, eqs. (fT9|) . establish a linear relation between the four left amplitudes A\ and the 
four right amplitudes 

A\ = M,,A] (20) 

where Mjj is the 4x4 "Matching Matrix" , not to be confused with the scattering matrix to 
be introduced later. 



B. Subsonic region 

Let us first find the roots of the dispersion relation f[T^ in the subsonic (right) region. 

This is graphically displayed in Fig. 8, where the sohd hne corresponds to the positive 

branch, the dashed line to the negative branch. 

So, for w > Q two real solutions exist, both belong to the positive norm branch, namely 
with positive group velocity (propagation towards the right, along the flow) and /c^ with 

negative group velocity (propagation towards the left, against the flow). These roots admit 
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Figure 9: Dispersion relation in the supersonic (left) region. 



a perturbative expansion for small w as 



hi, 



Vo - Cr 
UJ 

Vq + Cr 



1 + 777 ^ + Oiz^, 



(21) 



where the dimensionless expansion parameter used is Zr = The order terms are just the 
hydrodynamical modes discussed before. The other two solutions of the dispersion relation 
are a couple of complex conjugate roots k^. The one with positive imaginary part (/c!j_) 
represents a decaying mode as a; — )■ +oo. The other, k'^_, with negative imaginary part, a 
growing mode. These two complex roots do not exist in the hydrodynamical approximation. 
Their presence is a consequence of the nonlinearity of the dispersion relation which at /c ^ ^ 
becomes quadratic (single particle excitation). /c!j_ have a small w expansion 



cl - vl 



4(c2 - v^y 



2i 



cl-vl 



Having the four roots, we can write down the solution for the modes equations 
subsonic (right) region as 



(22) 
in the 



(23) 



and similarly for y^J^ (replace — )■ 



C. Supersonic region 



Let's move to the supersonic (x < 0, left) region. The dispersion relation (see Fig. 9) 
shows a completely different pattern, as compared to the subsonic case. 
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One sees that there exists a maximum frequency Wmax above which the description mimics 
exactly what happens in the subsonic case: two oscillatory modes /c^ and belonging both to 
the positive branch and propagating to the right and to the left respectively and two complex 
conjugate roots k^^. Note however that, although we are in the supersonic region, i.e. in 
the white hole, the kl mode is able to propagate to the left, i.e. upstream, notwithstanding 
the BEG is flowing supersonically in the right direction. The group velocity of this mode is 
clearly bigger than vq and also see from Fig. 9 that the corresponding root lives well outside 
the linear (hydrodynamical) regime. Wmax is the value of w for which ^\k^^^ = where 



kn 



1 




1/2 



(24) 



So for w > Wmax the solution of the mode equations ([9]) can be given in a way completely 
analogous to what was done in the subsonic case, namely eq. fl23|) . One can then impose 
the jump conditions fITI?]) on the left and right solutions across the discontinuity. We shall 
not discuss this case further since it just leads to a simple scattering process for the mode 
as would happen in a sub-sub configuration. 

From now on we shall restrict our discussion on modes with w < Wmax which will take us 
to a much richer physical outcome. Looking at Fig. 9 one sees that there are now four real 
roots A;^, A;[, fcg, k\ each corresponding to an oscillatory mode. The first two are the familiar 
hydrodynamical-like modes 



kl 



kl 



vo - Ci 
u 

Vo + Ci 



CiW^^f 



+ 0{zf) 
0{zf) 



(25) 



where zi = —. But while the k^ mode has, as before, positive group velocity and belongs to 
the positive norm branch, the k^ one, unlike what happens in the subsonic case, has positive 
group velocity and furthermore it belongs to the negative norm branch. The corresponding 
excitations are negative energy {w < 0) phonons dragged by the fluid to the right, towards 
the horizon, they are no longer able to propagate against the flow. The other two roots k'^, 
k[ are not hydrodynamical modes 



k 



3(4) 



^0 - cf 



\^iA±^M^^oizf) 



2 _ r^Yi 



2v^ 



1 - 



8{vi - c}f 



(26) 

These two additional oscillatory roots are the analytical continuation of the decaying and 
growing complex roots {k±) present in the subsonic regime (see eqs. (122|) ). These two modes 
have both negative group velocity, i.e. notwithstanding the supersonic character of the flow, 
they are both able to propagate against the flow (i.e. to the left). Clearly their group 
velocity is bigger than vq {\vg\ > Vq). Finally note that /C3 belongs to the positive norm 
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branch while ^4 to the negative one. Given this, the general solution of the modes equation 
in the left (supersonic) region can be given as (for w < Wmax) 



-ILUt 



(27) 



A similar expansion holds for ip^ with — )■ while the amplitudes A* remain the same 
as stressed before. 



D. The matching matrix Mij 

The four jump conditions at a; = (eqs. f|T9|) ) can be written as a linear system 



Wi 







/a:\ 


Ai 




K 




= Wr 




A'. 




A\ 


\A{) 




\A'i) 



where the 4x4 matrices Wi and read respectively 



W, 



Multiplying both sides of eq. 
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(28) 



(29) 



(30) 



(31) 



where M = W, ^Wr is the matching matrix. 
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Figure 10: Mode (pt^ri^)- 

E. The "in" scattering basis 

We now proceed, starting from the modes discussed in the previous sections, to construct 
a complete and orthonormal basis, with respect to the scalar product eq. f ITU]) . to be used to 
expand the field operator as indicated in eq. (|8]). This can be done in two ways. Either we 
choose incoming modes to build up the basis, i.e. modes that in the remote past {t — )■ — oo) 
start their journey from left or right spatial infinity {x = — oo, x = +oo) and propagate 
towards the discontinuity at x = 0, or we choose outgoing modes that at t = +oo propagate 
towards x = ±oo. 

For w > Wmax we have two incoming modes (A;^, k^) that are scattered (transmitted and 
refiected) in two outgoing modes {k^, kl). As said this case is trivial nd we shall not discuss 
it. 

For w < Wmax there are three incoming modes (/c^, k^, k^) and three outgoing modes 

We begin with the construction of the "in" basis. We define the scattering mode as 

following: an initial left moving unit amplitude v mode coming at t = — oo from the right 
subsonic region (i.e. /}jje~™*+*'^^'^) is scattered by the discontinuity at x = and generates a 
refiected right moving k^^ mode with amplitude (i.e. Al^^jj^e-^^^+^K^^ in the right region, 
and two transmitted k'^, k\ modes, with amplitudes and A\ respectively, in the left region 
(i.e. /l^/}^e-«"'*+«'«3^ and A\D'-^e~™*^'^^^'^^) . The construction is not complete since we need to 
include in the right (subsonic) region the complex decaying mode /c!|_ (i.e. y4!j_L'!|_e~™*^*'^+^). 
The growing mode is not included as it blows up at x = +oo. The construction of the 0^"(w) 
basis vector is depicted in Fig. 10. 

We stress again that k^ and k^ propagate upstream in the supersonic region and that /c4 
belongs to the negative norm branch. It represents the so called "anomalous" transmission 
and describes, as said, negative energy phonons. The associated amplitudes of the modes 
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Figure 11: Mode 



can be evaluated using the matching matrix M by solving the following system 

/ \ / 1 \ 





V4/ 



M 



\ / 



(32) 



The leading order results for small w read (we have set for simplicity % = m = 1) 

Vq - Cr 



u — 



^3 — 



/I4 _ 



'~-'ur,vr ) 
Vo + Cr 

w\/cl - vl[c] - cj) 



cl - vl + i\ vl - c^i) = S: 



3l,vr 



ALvr 



-s. 



3l,vr 



(33) 



The notation introduced is quite simple: the second couple of indices describes the incoming 
channel (vr) while the first couple of indices labels the outgoing channel. We will not need 
the exlicit form of A^. The amplitudes of the propagating modes can be shown to satisfy 
the unitarity relation 

I^;i' + I4l'-I4l' = l • (34) 

The minus sign in front of the /I4 term in eq. 
has negative norm. 

The second basis vector, the scattering mode 
depicted in Fig. 11. 

It corresponds to a unit amplitude positive norm right moving u mode coming from the 
(left) supersonic region (i.e. x^^e"*""*"*"*^"^) which is reflected into a positive norm /cg mode 
with amplitude (i.e. AgDge"™*"'"*'^^^) and a negative norm mode with amplitude A^^ 



341) comes from the fact that the kA mode 



[w), is constructed following the process 



(i.e. AiDie-'"^ 



In addition there is a transmitted mode in the subsonic region 
with amplitude moving to the right (i.e. ^^/}^e~™*"'"*'^^^) and a decaying mode with 
amplitude A^ (i.e. A!, L);e-™*+*'=+^). The amplitudes for this process are found again using 
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Figure 12: Mode (plliw). 



the matching matrix M 
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The leading order in w results are 
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ALul 



(36) 



The unitarity relation reads now |v4^P + |^3p — l^iP = 1- 

The last scattering mode, completing the "in" basis, is (p^^^iw) constructed according to Fig. 
12. 

An incoming unit norm, negative norm, mode (i.e. jj^^^-^^+^'^i^"^ jg reflected in a posi- 
tive norm /cg mode and in a negative norm ^4 mode with amplitude and respectively, 
and transmitted in a /c^ mode with amplitude AJ^. Furthermore, there is as usual in the 
subsonic region the decaying mode with amplitude A^_^. These amplitudes are found 
using the matching matrix as follows 



/ 1 \ 



\AJ 



M 



( \ 

K 
K 
\ / 



(37) 
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Figure 13: Mode Cr 



yielding 
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Cr ^^0 - Cl 



ur,vl 



Cl Vq + Cr 

(Vj - Cfr/\V0 - Cr) 
^2CiW^cl -vl[Ci - Cr 

A\ ~ -4* ^ S.,1 



At ^ 

^3 — 



.2 _ ,,2 

^0 



S: 



3l,vl 



The unitarity condition now reads 



I A 



r |2 



I A 



« |2 



lA 



« |2 



-1 . 



(38) 



(39) 



The minus sign on the r.h.s. comes from the fact that the unit amphtude incoming mode 
has negative norm. 

We have now the three "in" basis vectors (0^"(w), 0^'/(w), 'Pldi''^)) can expand the field 
operator as follows 



du 



^ in] 



(40) 



Note in this expansion that the mode 0^"(w), being a negative frequency one (negative norm) 
is associated to a creation operator, while (pZi'^) 0m" ('^) being positive frequency modes 
are associated to annihilation operators as usual. 

One can alternatively proceed to the construction of the "out" basis, which consists of unit 
amplitude outgoing modes (i.e. modes propagating away from the discontinuity at x = 0). 
These three modes have wave vectors /cg, ^4 and k"!^. Starting from each of these one can 
construct a corresponding basis vector. 

We begin defining the scattering mode 0°"*(w) constructed according to Fig. 13. 

This scattering mode is a linear combination of initial fc^, k^ right moving modes coming 
from the supersonic region and a left moving k^ mode coming from the right region with 
amplitudes A^^, Al, A^ respectively, which together with the decaying mode (of amplitude 
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Figure 15: Mode (t)"^\w). 



A^) produce a final right moving /cj^ mode in the subsonic region with unit amphtude. The 
corresponding amphtude for the incoming modes can be evaluated using the matching matrix 
M with a procedure identical to the one used for the "in" basis. Since we shall not need 
these amplitudes we shall omit the derivation of their explicit form. 

The other two scattering modes, completing the "out" basis, (l)'^^{w) and (f)°^^{w), correspond 
to the process depicted in Fig. 14 and Fig. 15 respectively. 

Hence the field operator can be expanded alternatively in the "out" basis as 



doj 



+ air + d^'\w)^'s'\w) + dt{w)^t\^)\ ■ (41) 

As before, note that the negative frequency (j)2f^{w) scattering mode is associated to a 
creation operator. 

Now, being both the "in" and "out" basis complete, the "in" and "out" scattering modes 
are linearly related, namely 



'•^"(U;) = Sur,vl<l>°urM + ^S^.i^sf (^) + S,i,,i(l)lf (w) . (42) 



This leads to a nontrivial Bogoliubov transformation between the "in" and "out" creation 
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and annihilation operators 

a°f{w) = S3i,vrat{w) + S3i,uiat:}{w) + S3i,viaZ\w) , 

= Sii^yraZ{w) + Sii^uiaZ{w) + Sii^yi%i^{w) . (43) 

The 3x3 scattering matrix 5 whose elements are the Sij^rni coefficients of eqs. f H3|) relate 
the "in" and "out" basis. One can show that it satisfies S'^riS = r] = SrjS^ where the 3x3 
matrix rj is t] = diag{l, 1, —1). 

As clearly seen from eqs. f H3|) . the scattering matrix mixes creation and annihilation opera- 
tors as a consequence of the fact that the basis contains both positive and negative frequency 
states and they get mixed passing from one basis to the other (see eqs. ( 142|) ). Because of 
this, the Hilbet space if^, constructed out of the "in" vacuum |0,m) (i.e. the state such 
that a'^j{w)\0,in) = for every by the action of creation dlj\w) operators, and the 
corresponding Hout Hilbert are not unitarily related. In particular \0,in) ^ \0,out). This 
means that, if we prepare the system to be at t = — oo in the "in" vacuum state, |0,m), i.e. 
the state containing no incoming quanta, at late time (t = +oo) this state (remember we 
work in the Heisenberg picture) will in general contain outgoing quanta. 
The physical effect we are therefore describing is the spontaneous emission of Bogoliubov 
(quasi) particles out of the vacuum that shows up in our EEC with a white hole like flow. 
This proceeds by the conversion of vacuum fluctuations into real on shell particles. The 
number of the different outgoing particles created out of the vacuum can be computed using 
eqs. as 

(0,m|a°ft(w;)a°f (u;)|0,m) = \Sur,vi\'' , 
{0,zn\dlf{w)a'^\w)\0,zn) = \S3iM\\ 

{0,tn\dlf{w)dlf{w)\0,in) = |54z,..|' + \S^i,ui\^ . (44) 

Here n'^ is (for given w) the number of /c^ particles emitted in the right (subsonic) region and 
travelling towards x = +oo, and n'^ are respectively the number of k^^ and k{ particles 
emitted in the left (supersonic) region: they travel supersonically with \vg\ > Vq against the 
flow towards x = —oo. Note that the /C4 quanta correspond to negative energy particles, 
they are the so called "partners" . From unitarity it follows 

<■ + nl = nt , (45) 

i.e. the number of positive energy particles equals the number of the negative energy ones 
as required by energy conservation. It is the existence of these partners in the spectrum of 
the physical states , that allows particles creation to occur even in a stationary setting, like 
ours, without contradicting energy conservation. 
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Figure 16: White hole configuration. 
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From eqs. 



Figure 17: Black hole configuration, 
we can infer the leading (in w) behaviour of these quantities 

31 _ {Cr - Vo){vi - C^y/^Cr + Ci) 1 
2ci{Cr + Vo){Cr - Cl) W 
41 _ {Cr - Vo){Vo-cf Y^^{Cr + Ci) 1 
2ci{Cr + Vo){Cr - Cl) W 



n. 



n. 



+ const. 
+ const. 



(46) 



F. White hole versus black hole emission 

In this section we shall discuss in detail the results for the particles emission obtained in 
the previous section for our white hole configuration and compare them to the ones obtained 
for a black hole configuration. The two are simply related: reversing the direction of the 
BEC flow, one simply moves from a white hole configuration to a BH one as shown in Figs. 
16, 17. 

For the black hole case sound waves are trapped in the x < supersonic region and 
dragged by the flow to the left, towards a; = — oo. The plot of the dispersion relation for the 
two configurations 



Wo — Vok = ±\ c^k"^ + 



(47) 



is represented in Figs. 18-21 . 
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Figure 18: Supersonic dispersion relation (left region) for the white hole. 
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Figure 19: Subsonic dispersion relation (right region) for the white hole. 



Moving from the WH configuration to the BH one {vq — > —Vq) one sees from Figs. 22, 
23 that the ingoing and outgoing character of the modes are exchanged, since BH and WH 
are related by time reversal and this operation reverses the sign of the group velocity. 

They are connected by a /c — j- —k symmetry as expected, since the dispersion relation eq. 
(H7|) is invariant under the combined k — )■ —k , vq — ?■ —vq reflection. 

Looking at the expression for the roots ki{w) for the BH case, given in the Appendix, and 
performing the above symmetry operation, one has the following correspondence 

k^lwH < — > k^lsH , 

kilwH < — y k[\BH ■ (48) 



At the fundamental level of the solution of the basic dynamical equations 
lish the correspondence between a WH and a BH mode as following 11 1 



one can estab- 



6f^(t,x)^0r''(t,a^) = [0r(-t,a:)]* 



(49) 



Figure 20: Supersonic dispersion relation (left region) for the black hole. 
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Figure 21: Subsonic dispersion relation (right region) for the black hole. 
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Figure 22: 'in' and 'out' modes in a white hole. 
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Figure 23: 'in' and 'out' modes in a black hole. 
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since under time reversal (t — )• —t, vq — )■ —Vq) eqs. dHD get mapped in their complex 
conjugated. The correspondence between the indices i and j is the one shown in (H8|l . 
Hence the WH "in" scattering modes ^^(j)^^{w) constituting the "in" basis for the field 
operator expansion can be related to the "out" BH scattering modes ^^(j)°^'^{w) (i.e. the 
"out" basis for the BH case) as follows 







X, 


vD = 






t, X, 








X, 


vD = 






t, X, 








X, 


vD = 




Kfi 


t, X, 





(50) 

where Vq^ = —v^^. This implies the simple relation between the WH scattering matrix 
SwH and the corresponding BH one Sbh 

[SwhT = [Sbh]-' , (51) 

which can re rewritten as 

SwH = vlSnnf V (52) 

where use of the unitarity of the scattering matrix {S'^riS = rf) has been made. 
Having established the relation between the S matrices for WH and BH fiow, we can now 
start to discuss our WH results of the previous section. As we have seen (eqs. dSl)) we 
have production of every three kind of particles. We begin with the emission into the (right) 
subsonic region, which can be considered as the exterior of the WH. Here we have emission 
of (m, r) quanta according to the first of eqs. 

WH UT_\c; |2 _ v^o ^Lj /^o\ 



ci {v^^ 



where we have added a superscript "TV if" to emphasyze that the fiow corresponds to a WH. 
The relevant channel responsible for this emission is in Fig. 24. The associated channel for 
a BH, according to the correspondence of (148|) is given in Fig. 25 and has an amplitude (see 
Appendix) 

cBH _ fCr '^O^ + Q _ /c7 Ci — _ ^ -^yn 

"''^^ V Cl Cr - \l ClCr+ V^^ ^ ^ 

in accordance with flS^ . 

Now eq. (153|) reveals that the emission from a WH in the subsonic region has at small 
w a fiat spectrum. It has not the 1/w thermal behaviour that characterizes the Hawking 
like emission from a BH fiow. As seen from eq. fl5^ this WH /c^ emission is related to the 
emission of negative energy fc^ "partners" inside a BH. 

Hawking radiation from a BH is described by the channel in Fig. 26 and the corresponding 
S-matrix element (see Appendix) gives 
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Figure 24: Relevant channel for the emission in the white hole exterior. 
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Figure 25: Associated black hole channel to that in Fig. 24. 



which describes a (approximately) thermal emission at an effective temperature 



Being (see eq. fl52 



T 



CBH 



,BH2 



,BH 



+ Cr 



CWH 



,BH\ 



(56) 



one sees that in a WH the process associated to 



Hawking like thermal emission occurring in a BH is the production of k''^ negative energy 
"partners" according to the channel in Fig 27. 

Finally, the production of positive energy /cg quanta inside a WH, Fig. 28, with = 
l^sfvil"^ (see the second of eqs. (jS]) and (H6|l ) is related to the following BH process Fig. 29, 
namely the production of kl partners inside the BH described by 5*^% with S^i^i = {—S^^i)* 
according to eq. (15^ . 



Conclusions 

We have seen that although the scattering matrices for BECs BHs and WHs quasi particle 
production are simply algebraically related (5'^^;^ = S]^]^), the physical features of the 
corresponding emission are rather different in the two contexts. 

The emission by a BEC WH in the subsonic region (i.e. outside the 'antitrapping' region 
bounded by the sonic horizon) has a flat spectrum that contrasts the thermal character of 
the BH emission in the subsonic region, which is the analogue of the process discovered 
by Hawking in the gravitational setting. A similar process indeed exists for BEC WHs 



25 



super 



=0 sub 



BH 



Figure 26: Channel responsible for black hole emission. 

> 



^0 



super sub 
WH 

Figure 27: White hole channel associated to that in Fig. 26. 

but occurs inside the horizon. The close relation existing between the two processes is 
summarized in S^^^i = —S^Jf^. However the produced quasiparticles are characterized by 
quite different modes in the two hydrodynamical 'u' mode in the BH case and a 

dispersive '^4' mode in the WH one. 

This leads to striking different physical features of the corresponding emissions. In par- 
ticular the associated density correlations functions inside WHs show a peculiar growing 
checkerboard pattern Q which has no analogue in the BH case (unless one considers a 



nonvanishing momentum transverse to the flow [12|). The responsible of this instabil- 
ity is a zero frequency propagating outgoing mode which is absent in the BH radiation 
(lim^^o K = 0' lim^^o K ^ 0) as can be seen in Figs. 18, 20. 

What can we infer from our analysis with respect to gravitational WHs (if these objects 
indeed exist)? First one should say that particle production by gravitational WHs has 
not been a subject much investigated leading however to controversial results. Just after 



Hawking's discovery in '74 of BH radiation Wald 13| analyzed the WH context concluding 



that the emission becomes singular. On the other hand, Hawking in IJ] claims that BHs 



and WHs behave exactly in the same way, i.e. also WHs should emit thermal radiation. 
Other claims that WHs 'antievaporate' are also present in the literature (see for instance 
15|). One should say that the results, unlike the BH case, seem to depend heavily on the 
WH formation process. Not to say about the role played by the past singularity and the 
boundary condition on it. 

Unfortunately our work like others on BEG WHs pseudo particles radiation does not 
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Figure 28: Channel associate to the production of positive energy /cg quanta inside a white hole. 
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Figure 29: Black hole channel associated to that in Fig. 28. 



shed light on the gravitational counterpart . The reason for this is that, as we have seen, 
the thermal like emission of BEC WHs is associated to the modes (/cg /C4 which 
are highly dispersive. The associated quanta differ significantly from relativistic massless 
quanta which are characterized by a non dispersive linear relation between frequency w and 
wavevector k. The ^3^4 modes are not hydrodynamical modes. The situation is just the 
opposite of what happens in BEC BHs thermal emission whose modes /c^ , k}^^^ are 
hydrodynamical and obey linear dispersion relation. 

Hence WH BECs thermal emission is, unlike the BH case, highly sensible to the disper- 
sion. One can say that it feels the short distance features of the underlying quantum theory. 
In this sense the robustness of BH radiation 16|| does not extend to WHs. 

One should always remember that the so called analogy between gravity and some con- 
densed matter systems relies on the possibility of treating the condensed matter system 
under the hydrodynamical approximation where modes are non dispersive and obey a linear 
'relativistic' dispersion relation. In view of all this one should hardly consider BEC WHs 
pseudoparticle production as a reliable insight for the corresponding gravitational setting. 
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Appendix A: Black hole S-matrix coefficients 



In this appendix we recall the ki{w) modes for the black hole, in the right (r) subsonic 
and left (l) supersonic regions 
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and the S-marix coefficients: 
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for the left moving 'in' mode from the subsonic region in Fig. 30; 
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for the positive norm right moving 'in' mode, Fig. 31, and finally 
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Figure 30: Incoming mode from the right subsonic region. 
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Figure 31: Positive norm 'in' mode from the left supersonic region. 
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Figure 32: Negative norm 'in' mode from the left supersonic region, 
for the negative norm 'in' mode from the left supersonic region. 
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